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ABSTRACT 

We report on the results from our deep Chandra observation (120 ks) of the super- 
nova remnant (SNR) N49 and soft 7-ray repeater (SGR) 0526-66 in the Large Magel- 
lanic Cloud. We firmly establish the detection of an ejecta "bullet" beyond the south- 
western boundary of N49. The X-ray spectrum of the bullet is distinguished from that 
of the main SNR shell, showing significantly enhanced Si and S abundances. We also 
detect an ejecta feature in the eastern shell, which shows metal overabundances similar 
to those of the bullet. If N49 was produced by a core-collapse explosion of a massive 
star, the detected Si-rich ejecta may represent explosive O-burning or incomplete Si- 
burning products from deep interior of the supernova. On the other hand, the observed 
Si/S abundance ratio in the ejecta may favor Type la origin for N49. We refine the 
Sedov age of N49, rg e( j ~ 4800 yr, with the explosion energy Eq ~ 1.8 x 10 51 erg. 
Our blackbody (BB) + power law (PL) model for the quiescent X-ray emission from 
SGR 0526-66 indicates that the PL photon index (r ~ 2.5) is identical to that of 
PSR 1E1048. 1-5937, the well-known candidate transition object between anomalous 
X-ray pulsars and SGRs. Alternatively, the two-component BB model implies X-ray 
emission from a small (R ~ 1 km) hot spot(s) (kT ~ 1 keV) in addition to emission 
from the neutron star's cooler surface (R ~ 10 km, kT ~ 0.4 keV). There is a con- 
siderable discrepancy in the estimated column toward 0526-66 between BB+PL and 
BB+BB model fits. Discriminating these spectral models would be crucial to test the 
long-debated physical association between N49 and 0526-66. 
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1. INTRODUCTION 



The supernova remnant (SNR) N49 is the third brightest X-ray SNR in the Large Magel- 
lanic Cloud (LMC). The blast wave is i nteracting with d ense clumpy interstellar clouds on the 
SNR's eastern side (jVancura et al. 19 92: Ba nas et al. 19971). producing bright emission in optical , 
infrared, ultraviolet, and X-ray bands (IPark et al. 2003al: ISankrit et al. 2004 : IWilliams et al. 20061 : 
Rakowski et al. 20071 : iBilikova et al. 20071 : lOtsuka et al. 20101 ) , which revealed the shock structures 
in N49 down to sub-pc scales over wide ranges of the velocity (v ~ 10 2-3 km s _1 ), temperature 
(T ~ 10 2 ~ 7 K), and the interacting density (n ~ 1-10 3 cm -3 ). The blast wave in the west- 
ern side of the SNR appears to be propagating into a lower density medium than in the east 
(jDickel &: Milne 19981 ). Early studies estimated overall low metal abundances in N49, indicating 
that metal-ri ch ejecta might have been substantially intermixed w i th the interstellar m edium (ISM) 
of the LMC (banzieer Leibowitz 19851 : iRussell fc Dopita 199ol : iHughes et al. 19981 ) . 



The origin of N49 has been elusive (thermonuclear vs. core-collapse exp losion) . The inters tellar 

environment of N49, i ndicated by the p resence of nearby molecular clouds (IBanas et al. 19971) and 

young stellar clusters (jKlose et al. 20041 ) and its location within an OB association (jChu Kennicutt 198a ) 
generally support a core-collapse origin for N49 from a massive progenitor star. One model sug- 
gested that N49 is interacting with a dense shell of a Stromgren sphere created by a massive 
B-type progenitor (jShull et al. 19851 ). The positional c o incidence of a soft 7-ra y repeater (SGR) 
0526-66 within the boundary of N49 (jCline et al. 19821 : iRothschild et al. 1994! ) may also suggest 
a core-collapse origin from a massive star, if this SGR is the compact remnant of the SN that 
created N49. Howeve r , the physical ass o ciation between SNR N49 and SGR 0526-66 is uncertain 
(|Gaensler et al. 200ll : klose et al. 20041 : badenes et al. 20091 ). Based on the ASCA data, Hughes 
et al. (1998) found no evidence for surrounding ISM that had been modified by stellar winds from 
the massive progenitor. 

Nucleosynthesis studies of metal-rich ejecta i n SNRs using spatially-resolved X-ray spectroscopy 

are useful to cons train the progenitor's mass (e.g. jHughes et al. 20031 : IPark et al. 2003bl : lPark et al. 2004 : 
Park et al. 20071 ). Although bright X-ray features in N49 are apparently dominated by emission 
from the shock-cloud interaction, previous Chandra observations have shown marginal evidence 
for enhanced emiss ion lines and overabundant metal elements from a few sub-regions in the SNR 
(jPark et al. 2003al ). For instance, a promising ejecta candidate was a small protrusion beyond the 
SNR's southwestern boundary. This emission feature was suggested to be a can didate ejecta bullet , 
but metal overabundances were not conclusive because of poor photon statistics (jPark et al. 2003al ). 
Also, because N49 was detected significantly off-axis (by ~6f5) in the previous Chandra data (Ob- 
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sID 1041) used by Park et al. (2003a), the point spread function (PSF) was substantially distorted, 
and the true morphology of this feature (point-like vs. ex tended) was uncertai n. We note that 
there were three other on-axis Chandra observations of N49/*1(: Kulkarni et al. 2003 ) . These previous 



observations were not useful to study the ejecta bullet candidate because of either the use of a 
CCD subarray or a very short exposure. Thus, a firm detection of metal-rich ejecta in N49 remains 
elusive. 

Here we report on initial results from our new Chandra observation of N49. Our new Chandra 
data present ~4 times deeper exposure (for the entire N49) than the data used by Park et al. 
(2003a), with an on-axis pointing to N49. With these new data, we clearly detect some metal-rich 
ejecta features in N49. Combining the new data with archival data, we also perform a detailed 
spectral analysis of SGR 0526-66. In this paper, we report on the results from our spectral analysis 
of SGR 0526-66 and the detected metal-rich ejecta features in N49. Our studies of complex shock- 
cloud interaction regions of N49 and the timing analysis of SGR 0526-66 will be presented in a 
subsequent paper. We describe the observations and the data reduction in Section [2j The image 
and spectral analyses are presented in Sectional and a discussion in Section HI A summary will be 
presented in Section [5j 



2. OBSERVATIONS & DATA REDUCTION 

We observed N49 with the Advanced CCD Imaging Spectrometer (ACIS) on board Chandra 
X-Ray Observatory on 2009 July 18 - September 19 (ObsIDs 10123, 10806, 10807, and 10808) 
during AO10. We chose the ACIS-S3 chip to utilize the best sensitivity and energy resolution of 
the detector in the soft X-ray band. The pointing (a 20 oo = 05 h 25 m 58 s . 8, <5 2 ooo = -66° 05' OO'.'O) 
was roughly toward the geometrical center of N49. SGR 0526-66 and the ejecta bullet candidate 
are detected within ~30" and ~40" of the aim point, respectively. We used a 1/4 subarray of the 
ACIS-S3 to ensure low photon pileup (<; 5%) of SGR 0526-66, and to allow us to detect the 8 s 
pulsations from the SGR, while still obtaining full coverage of the entire SNR N49. We processed 
the raw event files following the standard data reduction methods using CIAO 4.2 and Chandra 
CALDB 4.3.0, which includes correction for the charge transfer inefficiency (CTI). We applied 
the standard data screening by status and grade {ASCA grades 02346). We performed this data 
reduction on individual observations. The overall light curve from the entire S3 chip for each ObsID 
showed no evidence of flaring particle background. Then, we merged four data sets into a single 
event file by re-projecting them onto ObsID 10807's tangential plane. After the data reduction, the 
total effective exposure is ~108 ks. 

We also used two data sets (ObsIDs 747 and 1957) available in the Chandra archival The 



ObsIDs 747, 1957, and 2515 with exposures of 40, 48, and 7 ks, respectively. 

2 In the Chandra archive, there are two other data sets (ObsIDs 1041 and 2515) that detected N49 and SGR 
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pointing of these archival data was toward the position of SGR 0526-66 (a 20 oo = 05 h 26 m 00 s . 7, 
^2000 = —66° 04' 35"0), and thus they provide a high resolution imaging of N49 as well as the 
SGR. Although these archival data cannot be used to study south-southwestern regions (including 
the ejecta bullet candidate) of N49 due to the use of 1/8 subarray of the ACIS-S3, they covered 
east-northeastern parts of the SNR with a decent exposure (a total of ~88 ks by combining ObsIDs 
747 and 1957). Thus, these data are useful to study eastern regions of N49 and SGR 0526-66. We 
re-processed these two archival data sets in the same way as we did for our new observation data. 
Combining all these observations, our data present the deepest coverage with the high resolution 
X-ray imaging spectroscopy for SNR N49 (~108 ks in the south-west regions and ~196 ks in the 
north-east regions) and SGR 0526-66 (~196 ks). Our new Chandra observations and archival data 
used in this work are summarized in Table [TJ 



3. ANALYSIS & RESULTS 



3.1. X-ray Images 



An X-ray color image of N49 shows complex asymmetric X-ray emission features ( Figure [HQ . 

These features have been partially revealed by previous Chandra data (jPark et al. 2003al ; iKulkarni et al. 20 03) 
Our new deep Chandra image reveals all those fine X-ray structures throughout the entire SNR 
with substantially improved photon statistics and resolution (Figure Hh) . The outer boundary of 
N49 appears as a thin shell radiating in soft X-rays (reddish in Figure [Th.) that represents the blast 
wave sweeping through the general ambient medium of the LMC. A notable exception is the can- 
didate ejecta bullet extending beyond the SNR shell in the southwest, which is distinctively blue 
in color (Figure [lb). Bright X-ray filaments in t he eastern parts o f the SNR are regions where the 
shock is interacting with dense clumpy clouds (jPark et al. 2003al . and references therein). These 
shock-cloud interaction regions show a wide range in X-ray colors at angular scales down to several 
arcseconds. The western parts of N49 are relatively faint and X-ray emission is emphasized in 
green to blue while showing some red filamentary features as well (Figure [T^,) . This general color 
variation across the SNR appears to be related to the temperature variation caused by the overall 
interstellar density distributio n, in which dense c louds are interacting with the shock mostly in the 
eastern half of the SNR (e.g., IPark et al. 2003ah . 



The spectrally-hard, blue protrusion extending beyond the main SNR shell in the southwestern 
boun dary was proposed to be a candidate ejecta bullet based on the previous low-resolution Chandra 
data (jPark et al. 2003a!). If it is indeed an ejec ta bullet, this feature is reminiscent of "shrapnel" 
found in the Vela SNR (jAschenbach et al. 19951 ). but nearly two orders of magnitude more distant. 
Our new Chandra data clearly resolve this feature. It is an extended feature with a head (~4" in 



0526-66. However, these data are not useful for this work because of a large off-axis angle for N49 (~6l5, ObsID 
1041) and a short exposure (~7 ks, ObsID 2515). 



- 5 - 



radius) and a tail that connects the head to the main SNR shell (the "Head" and "Tail" regions in 
Figure QJi, respectively). The blue color of the Head region is likely due to enhanced line emission 
in the hard X-ray band: e.g., the Si (Hea + Lya) line equivalent width (EW) marjf] shows a strong 
enhancement in the Head region (Figure [2k.) , which suggests a metal-rich ejecta for the origin of 
this feature. In fact, our spectral analysis of the Head region shows highly overabundant Si and S 
(see Section [3.2.2j) . The Tail region appears to be composed of a central hard (blue) emission region 
surrounded by a soft (red) , conical region (Figure QJ) . It is likely the turbulent region behind the 
bullet, in which the metal-rich ejecta and shocked ISM are intermixed, possibly enclosed by the 
sides of a bow-shock produced by the bullet's supersonic motion. The bow shock, however, is not 
prominent ahead of the bullet, possibly due to unfavorable viewing geometry. The axis bisecting 
the opening angle (~87°) is roughly pointing back to the geometric center of the SNR, but not 
exactly toward SGR 0526-66. 

We noticed that there is a small region in the eastern part of the SNR (the "East" region in 
Figure Hh), which shows similar characteristics to those of the Head region: i.e., this region stands 
out with a distinctively blue color in the bright eastern half of the SNR, and it is coincident with 
a strongly enhanced Si line EW (Figure [2h.) . The East region corresponds to an inter-cloud region 
(Figure [2b) in wh ich a possibility of metal overabundance was suggested by our previous work 



(jPark et al. 2003al ). Our spectral analysis of this region indeed reveals its metal-rich ejecta nature 



(see Section [3.2.2[) . 



3.2. X-ray Spectra 

Based on results from our image analysis, we extracted spectra from several regions that 
appear to characteristically represent distinctive features: i.e., the swept-up interstellar medium 
(the "Shell" and "SGR_BG" regions), shocked metal-rich ejecta (the "Head" and "East" regions), 
and the connecting region between the ejecta bullet and the main SNR shell (the "Tail" region). 
These regions are marked in Figures QJ, &[lb. We also extracted the spectrum from the quiescent 
X-ray counterpart of SGR 0526-66 (the "SGR" region in Figure [lb). We present the spectral 
analysis of these swept-up medium, ejecta features, and 0526-66 in Sections 13.2. 1\ [3T2T21 and l3.2.3l 
respectively. 



3 We created the Si line EW map following the method described in Park et al. (2003a). The 1.75-2.1 keV band 
was used to extract the Si line image. The upper and lower energy continuum images were created using the 2.14-2.28 
keV and 1.63-1.71 keV bands, respectively. 



- 6 - 



3.2.1. N49: Swept- Up Interstellar Medium 
The brightest X-ray emission in the eastern parts of N49 originates from t he shock interaction 



with clumpy clouds with a density gradient spanning a few orders of magnitudes (jVancura et al. 19921 ) 
X-ray emission from such a clumpy medium would represent a complex mixture of gas temperatures 
and ionization states. Such a plasma condition is likely localized in strong shock-cloud interaction 
regions. While the western regions do not show evidence of shock interaction with dense clumpy 
clouds, they also show some significant X-ray color variation (Figure . For those regions, it is 
not straightforward to characterize the X-ray emission spectrum of the general ambient medium 
shocked by the blast wave. Thus, to study the nature of the general swept-up ISM and the SNR's 
overall dynamics, we extracted the X-ray spectrum from a region in the outermost rim of the SNR 
in the southern boundary (the "Shell" region in Figure [Th.) in which no complex spectral and/or 
spatial substructures are seen. We also used a small annular region surrounding SGR 0526-66 (the 
"SGR_BG" region in Figure [lb) to characterize background emission for SGR 0526-66. Since the 
SGRJ3G region is in the transition area between the bright eastern and faint western parts of the 
SNR, it may also effectively show an average plasma condition of the complex shock structures 
across the SNR. 

For our spectral analysis, we re-binned the observed spectra to contain a minimum of 20 
counts per energy channel. We fit the X-ray spectrum with a non-equilibrium ionization (NEI) 
state pla ne-shock model (vpsho ck model in conjunction wit h the NEI version 2.0 in the XSPEC 



software. iBorkqwski et al. 200ll ) that is based on ATOMDB (jSmith et al. 200 ll ). We used an aug 



mented versio^ of this atomic database to inclnde inner-shell processes lines from Li-hie 
ions) and updates of the Fe L-shell lines, whose effects are substantial in under-ionized plasma, 
but are not incorporated in the standard XSP EC NEI version 2.0. W e fixed the Galactic col- 



umn at iVH,Gai = 6 x 10 cm toward N49 (jDickey fe Lockman 19901 ). We fit the foreground 



column by the LMC (-/Vh,lmc) assuming the LMC abundances available in the literature (He = 
0.89, C = 0.30, N = 0.12, O = 0.26, Ne = 0.33, Na = 0.30, Mg = 0.32, Al = 0.30, Si = 0.30, 
S = 0.31, CI = 0.31, Ar = 0.54, Ca = 0. 34, Cr = 0.61, Fe = 0.36, Co = 0.30, and Ni = 0.62, 



Solar (jAnders fe Grevesse 19891 ). 



Russell fe Dopita 19921; iHughes et al. 19981 ). Hereafter, elemental abundances are with respect to 



The southern boundary of N49 was not covered by ObsIDs 747 and 1957 because of the use of 
a 1/8 subarray. Thus, we used our new Chandra data to extract the spectrum from the Shell region. 
The background was subtracted using the spectrum extracted from two circular source-free regions 
(with a radius of 8") beyond the south-western boundary of the SNR. The Shell region spectrum 
contains significant photon statistics of ~4600 counts, and can be well fitted by the NEI plane-shock 
model (x 2 / n = 49.1/68 with the electron temperature kT = 0.57 keV and the ionization timescale 
n e t = 6.35 x 10 11 cur 3 s, Figure [3^). The foreground column by the LMC is estimated to be 



4 This augmented NEI model has been provided by K. Borkowski. A relevant discussion on this modeling issue 
can be found in Badenes et al. (2006). 
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-^h,lmc ~ 0.9 x 10 21 cm 2 . The best-fit metal abundances are generally consistent with the LMC 
values. Results from this spectral model fit are summarized in Tables [2] & [3j 

The SGFLBG region was observed by ObsIDs 747 and 1957 as well as our new Chandra 
observations, and thus, we used all of those data to extract the spectrum from this region. The 
total photon statistics for this region is ~13600 counts. The background was subtracted using the 
spectrum extracted from a circular source-free region (with a radius of 8") beyond the northern 
boundary of the SNR. We initially fit these three data sets simultaneously with all fitted parameters 
tied with each other. Although the overall fit may be acceptable (x 2 / n = 1-28), there appears to 
be a small systematic offset in the normalization at E <^ 1 keV between data taken in 2000/2001 
(ObsIDs 747 and 1957) and in 2009 (ObsIDs 10123+10806+10807+10808). Considering the fact 
that the time separation between 2000/2001 and 2009 observations is substantial and that the effect 
appears to be emphasized in the soft band {E keV), this small discrepancy is likely related to 
the calibration effect by the time-dependent quantum efficiency degradation in the ACIS. Since the 
quantum efficiency degradation affects the spectral modeling as if there is an "extra" foreground 
absorption, we attempted the model fit with -/Vh,lmc untied among three spectra. The best-fit 
-^VhjLMC values are the same between 2000 and 2001 data, while it is ~25% larger for the 2009 
data. Thus, we repeated the spectral model fit with A^h.lmc tied only between 2000 and 2001 
data, and with that for 2009 data allowed to vary freely. The fit improved somewhat (x 2 / 11 = 
1.19). The best-fit columns are iV H ,LMC = l-39±g;^ x 10 21 cm" 2 and 1.76±g;|| x 10 21 cm~ 2 for the 
2000/2001 and 2009 data, respectively. While showing a marginal calibration effect, these values 
are consistent within statistical uncertainties (hereafter, errors are with a 90% confidence level). In 
the following discussion, we assume the average value of -/Vh,lmc = 1-58^q' 44 x 10 21 cm -2 . The 
best-fit spectral parameters for the SGR_BG region are kT ~ 0.56 keV and n e t ~ 9.65 x 10 11 cm -3 
s, all of which are consistent with those for the Shell region within uncertainties. The best-fit metal 
abundances are also generally consistent with those from the Shell region. These model parameters 
are summarized in Tables [2] & [3l 

3.2.2. N49: Metal-Rich Ejecta 

Our new Chandra data clearly resolve the extended nature of the ejecta fragment bullet (the 
Head region) and the trailing tail-like hot gas (the Tail region, Figure [It). Because the bullet region 
was not covered by the 2000/2001 data, we used only 2009 data for the spectral analysis of this 
feature. We extracted ~1800 counts from the Head region. In contrast to those for the swept-up 
medium, the X-ray spectrum of the Head region shows remarkably enhanced emission lines from 
highly ionized Si and S (Figure 0ti), which is consistent with the corresponding bright Si line EW 
image (Figure [2^i) . We fitted this spectrum with an NEI plane-shock model following the method 
described in Section 13.2. li The background was subtracted using the same spectrum as we used 
for the Shell region. We initially fixed metal abundances at the best-fit values for the Shell region 
(Table [3]), which resulted in a statistically unacceptable fit (x 2 / n = 3.3) primarily because of the 
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strongly enhanced Si and S lines. Varying Si and S abundances provides an acceptable fit (x 2 / n 
= 67.0/59, Figure Hi). The best-fit LMC column (iV H ,LMC ~ 1-3 x 10 21 cm- 2 ) is consistent with 
those estimated from the Shell and SGR_BG regions. The electron temperature (kT ~ 1 keV) and 
the ionization timescale (n e t > 2 x 10 12 cm -3 s) are significantly higher than those from swept-up 
medium. The Si (~2.3) and S (~3.2) abundances are an order of magnitude larger than the LMC 
values, firmly establishing the metal-rich ejecta nature of the bullet. 

Varying other elemental abundances (O, Ne, Mg, and Fe) somewhat improves the fit (x 2 /n 
= 49.3/55, but the corresponding F-probability is marginal [~0.002]), and the best-fit electron 
temperature is very high (kT ~ 4 keV). In addition to the overabundant Si and S, the best-fit 
abundance appears to be moderately enhanced for O (~0.7), Ne (~2.1), and Mg (~1.0). However, 
these abundances are not well constrained (with a factor of ^2 uncertainties), and thus the over- 
abundance for these elements is not convincing based on the current data. We conclude that the 
main differences in the fitted parameters resulting from varying metal abundances other than Si 
and S (i.e., the high kT ~ 4 keV and possibly enhanced abundances for O, Ne, and Mg) are not 
compelling. Otherwise, it does not affect our discussion of the ejecta features (see Section B~T]) . 
Thus, we hereafter discuss the nature of the bullet based on our spectral model fit with only Si and 
S abundances varied. The results from the Head region model fit are summarized in Tables [2] & [3j 

The Tail region is faint, and we extracted ~1500 counts in this region. We fit this spectrum 
with an NEI plane shock model x 2 / n = 44.9/44, Figure Hh). The Tail region is fitted by the X- 
ray spectrum from a hot gas (kT ~ 2 keV) in a significantly under-ionized state (n e t ~7 x 10 10 
cm -3 s). The best-fit elemental abundances are consistent with those estimated from the Shell and 
SGR_BG regions. We note that the Tail region shows substructures of the spectrally-hard (blue) 
interior surrounded by soft (red) outer layer (Figure [lb). It suggests that the Tail region may be 
a mixture of metal-rich ejecta and shocked swept-up medium. Thus, we attempted an alternative 
spectral model, assuming two shock components, a soft component for the swept-up ISM and a 
hard component for the metal-rich ejecta. The spectral parameters, except for the normalization, 
and metal abundances for the ejecta component were fixed at the best-fit values estimated from the 
Head region. The metal abundances for the soft component were fixed at the values estimated from 
the Shell region. The electron temperature, ionization timescale, and the normalization parameters 
for the soft component were varied freely. The observed Tail spectrum can be equally fitted by 
this two-component model (x 2 / n = 52.4/49). In this model fit, the metal-rich ejecta component 
contributes ~15% of the total observed X-ray flux in the 0.5-10 keV band. These one- and two- 
component models are not statistically distinguished with the current data. Results from these 
Tail region spectral model fits are summarized in Tables [2] & [3j 

The East region shows a distinctively blue color (Figure QJ,) with a strong enhancement in 
the Si line EW image (Figure [2^t) , which resembles the spectral properties of the H ead region of 
the bullet. This region corresponds to a hot inter-cloud region ()Park et al. 2003a! ) . This region 
may also be a candidate ejecta feature, probably isolated from the complex emission from shocked 
dense clumpy clouds (Figure [2b). Since this region was covered by both the 2000/2001 and 2009 
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data, we used all of these data to extract the X-ray spectrum. We extracted the X-ray spectrum 
from a small circular region with a radius of 2", in which the Chandra data allowed us total photon 
statistics of ~3600 counts. The background was subtracted using the spectrum extracted from two 
circular source-free regions (with a radius of 8") beyond the eastern boundary of the SNR. We fit 
the observed spectrum with the NEI plane shock model. Because enhanced Si-K and S-K lines are 
evident, we varied Si and S abundances while fixing other abundances at average values estimated 
from the Shell and SGR_BG regions. The fitted parameters are -/Vh,lmc ~ 0.8 x 10 21 cm" 2 , kT ~ 
1.05 keV, and n e t > 52.4 x 10 11 cm~ 3 s. The estimated Si and S abundances are high (Si ~ 1.8 and 
S ~ 1.4). While the fit may be acceptable (x 2 / n = 173.4/134), there appear to be some systematic 
residuals at E ^ 1 keV. Based on the results from our spectral analysis of the ejecta (the Head 
region) and the swept-up ISM (the Shell and SGFLBG regions) regions, we consider that the East 
region emission is a superposition of the spectrally-hard metal-rich ejecta and the projected soft 
swept-up ISM components. To account for the contribution of soft X-ray emission from the shocked 
ISM in the observed spectrum of the East region, we added an NEI model component for which we 
fixed the electron temperature, ionization timescale, and metal abundances at the average values 
estimated from the Shell and SGR_BG regions. Only the normalization parameter was varied for 
this component. This model fit significantly improves (x 2 / n = 144.4/133) without affecting the 
best-fit parameters for the hard, ejecta component. The only noticeable change is a larger value for 
the best-fit -/Vh,lmc ~ 3 x 10 21 cm~ 2 . This larger estimate for JVh,lmc appears to be consistent 
with the presence of dense molecular clouds toward eastern parts of the SNR. The assumed soft 
component from the swept-up ISM contributes ~20% of the total observed X-ray flux in the 0.5-5 
keV band. The results from this fit are summarized in Table [2] & [3l 



3.2.3. SGR 0526-66 

The X-ray spectrum of SGR 0526-66 is extracted from a circular region of 2" radius (Figured^). 
The background spectrum is characterized by emission from an annular region (2" in thickness) 
around the source (the SGR_BG region, Figure [lb) to average the complex diffuse emission from 
SNR N49 projected against the SGR. Since archival Chandra data (ObsIDs 747 and 1957) were 
pointed at the SGR with a decent exposure (Table [I]), we used these archival data as well as our 
new data for the spectral analysis of SGR 0526-66. We fit these three spectra simultaneously. Just 
like our spectral analysis of N49, we assumed that the SGR spectrum is absorbed by both Galactic 
(fixed at -/V"H,Gai = 6 x 10 20 cm -2 ) and LMC columns (A^h.lmc)- 

Initially, we fit the SGR spectrum with a single power law (PL) in which we tied all fitted 
parameters between the three sets of the spectrum. The best-fit model from this fit is not ac- 
ceptable (x 2 / n ~ 1.8). The best-fit model shows significant residuals at E J> 0.8 keV, because it 
underestimates the overall flux level for the 2000/2001 data while it overestimates the flux level 
for the 2009 data. Thus, we untied the normalization parameter between the 2000/2001 and 2009 
data, and then repeated the model fit. The fit improved (x 2 / n ~ 1-5) with a ~15% lower nor- 
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malization for the 2009 data than that for the 2000/2001 data. Although statistically improved 
by removing the normalization offset between the 2000/2001 and 2009 data, we do not accept this 
fit because systematic residuals are evident over the entire bandpass with a relatively high x 2 / n 
~ 1.5. Alternatively, we attempted a single blackbody (BB) model fit. The best-fit models are 
statistically unacceptable (x 2 / n ~ 2.2 and 2.5 with the untied and tied normalization parameter 
between the 2000/2001 and 2009 data, respectively). 

Kulkarni et al. (2003) showed that a PL+BB model was needed to adequately describe the 
quiescent X-ray emission from SGR 0526-66. On the other hand, Park et al. (2003a) argued 
that there could have been a considerable contamination from the surrounding soft thermal X-ray 
emission in the SGR spectrum presented by Kulkarni et al. (2003), and that the need for a BB 
component was not conclusive. Our SGR spectrum with ~4 times higher photon statistics than 
that used by Park et al. (2003a) indicates that a PL model cannot describe the observed SGR 
spectrum, supporting the conclusion by Kulkarni et al. (2003). Thus, we fit the SGR spectrum 
with a PL+BB model. We repeated the general steps as we did with a single PL model fits: i.e., we 
first tied all fitted parameters between the 2000/2001 and 2009 data sets, and then untied them one 
at a time (i.e., BB temperature [&Tbb], PL photon index [r], and normalizations for the BB and PL 
components). The best-fit values for /cTbb and T are consistent within uncertainties between the 
2000/2001 and 2009 data. The normalization for the BB and/or PL components are fitted to be 
~10-30% lower for the 2009 data than they are for the 2000/2001 data. The fit improvement due to 
this normalization change is statistically significant, and the fit is acceptable (x 2 /n = 577.4/582). 
The normalization change appears to be more effective for the PL component than for the BB, 
probably because the PL component covers a broader range of the observed spectrum than the 
BB component (The PL component contributes ~70% of the total observed flux in the 0.5-10 keV 
band). 

While the X-ray s pectra of SGRs an d anomalous X-ray pulsars (AXPs) are usually fitted with 



a BB+PL model fe. g.. iM ereghetti 2007) , doub le BB models (BB+BB) have also been suggested to 



fit some AXPs (e.g., lHalpern fc Gotthelf 20051 ). A recent work showed that a BB+BB model fit can 



successfully reproduce the observed Chandra spectrum (taken in 2000 and 2001) of SGR 0526-66 



(jNakagawa et al. 20091 ). Thus, we attempted to fit our SGR spectrum with a BB+BB model. We 
tied BB temperatures between 2000/2001 and 2009 data since we found no evidence of change in the 
temperatures. Just like our BB+PL modeling, we varied normalizations freely between 2000/2001 
and 2009. The fit is statistically good (x 2 / n = 616.4/582). Although statistical uncertainties are 
relatively large, a ~15% decrease in the best-fit value of the normalization parameters are suggested 
for both BB components, which is similar to the results from the BB+PL model fit. 

Both of the BB+PL and BB+BB model fits suggest that the overall X-ray flux of SGR 0526-66 
in the 0.5-10 keV band has been reduced by ~15% in 2009 compared with that in 2000/2001. We 
note that we used a 1/4 subarray of the ACIS-S3 in our 2009 observations, while a 1/8 subarray 
was used in 2000 and 2001. Based on the event grade distribution, we found that the photon pileup 
effect is similarly small in both data sets (~2% in 2000 and 2001 vs. ~4% in 2009). Nonetheless, 



- 11 - 



assuming the X-ray flux observed in 2000/2001, we estimated the pileup effect on the X-ray flux 
measurements in 2009 data using the Portable Interactive Multi-Mission Simulator (PIMMS). Our 
PIMMS simulations predict that the pileup effect in 2009 data would reduce the 0.5-10 keV band 
X-ray flux in 2009 by ^5% compared with t hat measur ed in 2000/2001. We also tested the pileup 
effect by applying an ACIS pileup model ([Davis 200ll ) for our spectral model fits. This model 
indicated that the 0.5-10 keV flux estimates with the 2009 data are probably affected by <^5% due 
to the use of 1/4-subarray. Thus, although the flux change between 2000/2001 and 2009 is partially 
caused by pileup effect, the pileup is not responsible for the entire flux change. The apparent flux 
change is unlikely due to uncertainties in the detector calibration for the time-dependent quantum 
efficiency degradation, because there is no evidence for discrepancy in the fitted normalization 
parameters that is emphasized in the soft band (E 1 keV) between 2000/2001 and 2009 data. 
Thus, while follow-up Chandra observations are required to draw a firm conclusion on the nature 
of the long-term X-ray light curve of 0526-66, we tentatively conclude that the discrepancy in the 
observed X-ray flux of SGR 0526-66 between 2000/2001 and 2009 is probably real rather than an 
artifact due to the photon pile- up and/or an inaccurate calibration of the time-dependent quantum 
efficiency degradation between two epochs. The SGR spectrum and the best-fit BB+PL model are 
presented in Figure EH. Results from the best-fit BB+PL and BB+BB model fits are summarized 
in Table H] and EJ respectively. 



4. DISCUSSION 
4.1. N49 

Based on the volume emission measure (EM) estimated from the spectral fit of the Shell 
region, we calculate the post-shock electron density (n e ) behind the blast wave. We assumed an 
emission volume of V ~ 9.4 x 10 55 cm 3 for a 3" x 18" region (corresponding to ~0.7 pc x 4.4 pc, 
hereafter we assume d = 50 kpc for the LMC) with a ~1 pc path-length along the line of sight. 
We also assumed n e = 1.2 nn for a mean charge state with normal composition. The best-fit EM 
(= 6.75 x 10 57 cm~ 3 ) implies n e ~ 9.3/~5 cm -3 and n# ~ 7.7/~2 cm~ 3 , where / is the volume 
filling factor of the X-ray emitting gas. The pre-shock H density is then no ~ 1.9/ _ 5 cm -3 for a 
strong adiabatic shock where nu = 4no- Assuming an adiabatic shock in electron- ion temperature 
equipartition, the gas temperature is related with the shock velocity (v) as T = 3fhv 2 /16k where 
k is the Boltzmann constant and fh ~ 0.6 m p is the mean molecular weight with the proton mass 
m p = 1.67 x 10~ 24 g. For the gas temperature of kT = 0.57 keV, the shock velocity v ~ 700 km 
s _1 is estimated. Thus, for the radius of ~8.5 pc (assuming the angular distance of ~35" between 
the Shell region and the geometric center of N49), we estimate the Sedov age of the SNR, rg c d ~ 
4800 yr. The estimated rg e d and no imply an explosion energy of Eq ~ 1. 8f~2 x 10 51 erg. The 



estimated Sedov age is ~30% lower than our previous est imate (^6600 vr. iPark et al. 2003al ) for 



which we assumed no = 0.9 cm 3 based on optical data (IVancura et al. 19921 ) and the canonical 



3) 



value of the SN explosion energy (Eq = 1 x 10 erg). Also, the previous estimate was based on 
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a more complex spectral modeling of multi-temperature plasma in the eastern regions of the SNR 
for which an accurate estimate of the distance from the SNR center was not feasible. Thus, we 
conclude that our new estimate of the SNR age is more reliable, while being in plausible agreement 
with the previous estimate by Park et al. (2003a). 

The presence of the ejecta bullet beyond the southwestern shell of N49 is conclusively estab- 
lished by our new Chandra observation. The on-axis ACIS image clearly resolves the bullet into a 
head and a tail extending back to the main shell of the SNR. The head region itself is an extended 
feature (~4" in radius) with an enhanced intensity toward the inferred direction of motion (Fig- 
ure [lb) . The X-ray color of the head is distinctively blue in contrast to reddish color of the overall 
SNR shell. The foreground column (JVh,LMC ~ 1-3 x 10 21 cm -2 ) toward the head is consistent 
with those to N49's main shell and the bullet's tail region, which supports that the blue color is 
intrinsic for the head rather than being caused by a significantly larger foreground absorption (in 
which case the head would probably be a background source). In fact, the observed X-ray spectrum 
(Figure Hti) shows that the hardness of the head is due to highly enhanced line emission from He- 
and H-like Si and He-like S ions. The estimated abundances for Si (~2.3) and S (~3.2) are an 
order of magnitude larger than the LMC values. We note that Park et al. (2003a) suggested a 
possibility of the enhanced abundance for O in addition to Si and S in the bullet. However, the 
overabundance for O is not confirmed by our new data. Nonetheless, the highly overabundant Si 
and S firmly establish the identification of the bullet as metal-rich stellar fragment ejected from 
the deep interior of the progenitor star. 

The East region shows nearly identical spectral characteristics to the Head region: i.e., a dis- 
tinctively blue color compared with surrounding regions and strongly enhanced Si and S lines. Just 
like the Head region, significantly enhanced Si and S abundances are estimated, but overabun- 
dances for other metal species are not required to describe the observed spectrum. The foreground 
absorption for the East region appears to be larger than other regions of N49, and it is probably 
cause d by nearby dense molecular clouds whi ch are interacting with the SNR in the eastern regions 



(e.g., iBanas et al. 19971 : lOtsuka et al. 2010i ). If the excess column (AiY H ,LMC ~ 2 x 10 21 cm -2 
compared with the average -/Vh,lmc ~ 1-2 x 10 21 cm -2 estimated for other regions of N49, Table[2]) 
originates from these nearby molecular clouds, the corresponding average cloud density of n# ~ 90 
cm -3 is implied for the overall cloud size of R ~7 pc (the size estimated by Banas et al. 1997). This 
average density is in plausible agreement with the pre-shock density range of the clumpy clouds (no 
~ 20-940 cm" 3 , Vancura et al. 1992) with which N49 is interacting, while it is significantly larger 
than the value we estimated for the Shell region (no ~ 1.9 cm -3 ) where the shock is propagating 
into the low density ambient medium. 

Based on the Si and S abundances and volume emission measures of these metal-rich ejecta 
features, we estimate the Si to S ejecta mass ratio. The observed spectrum of the Head and East 
regions shows that He- and H-like ionization states dominate for Si, while S ions are primarily in 
He-like state. Thus, for simplicity, we assumed a "pure" ejecta case with electron to ion density 
ratios of n e> si ~ 12.5 ns; and n ej s ~ 14 ns- For dominant isotopes of 28 Si and 32 S, the measured 
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Si and S abundances imply the ejecta mass ratio Mgi/Ms ~ 1.2 and 1.6 (Vjgi/Vs) 5 , where Vsi and 
Vs are the volume of Si- and S-rich hot ejecta gas, respectively. Assuming Vsi ~ Vs, the estimated 
mass ratios are Mgi/Mg ~1.2 for the Head region and ~1.6 for the East region. The mixture of 
H in metal-rich ejecta may not significantly affect our mass ratio estimates as long as the ISM 
mixture rate is similar between the Si and S ejecta. We compare these Mgi/Ms with standard 
SN nucleosynthesis models. Our estimated Mg;/Afg appears to be smaller than those for core- 
colla pse SN models in which Ms\/Ms typically r anges ^2-4 depending on the progenitor's mass 
(e.g., iNomoto et al. 1997al ; iRauscher et al. 20021 ; iLimongi Sz Chieffi 20031 ). The Si to S mass ratio 
for Type la SN models (M Si /M s ~ 1.5 - 1.8, e.g., Nomoto et al. [1997b]; Iwamoto et al. [1999]) are 
generally smaller than core-collapse cases, which are closer to our estimates for N49. The lack of 
evidence for O-rich ejecta in the Head and East regions is generally suggestive of a Type la origin 
as well. If we take ~0.3 (Table [3j) as an upper limit for the O abundance of the SN nucleosynthesis 
in N49, Mo/Msi < 1.4 can be inferred for the Head region. This limit appears to be consistent 
with Type la models (Mo/Mgi <^ 1, e.g., Nomoto et al. [1997b]; Iwamoto et al. [1999]), while being 
smaller than those for core-collapse models (Mq,/Mq,\ > 2, e.g., Nomoto et al. [1997a]; Limongi & 
Chieffi [2003]). 

On the other hand, the Head and East regions do not show evidence for overabundant Fe, 
which is usually considered as an iconic feature for Type la SNRs. The lack of Fe-rich ejecta is 
problematic for a Type la interpretation, and it would rather support a core-collapse origin. Also, 
the suggested Type la origin for N49 is inconsistent with its interstellar environment with recent 
star- formi ng regions and nearby molecular clouds, which rather suggest s a massive progenito r for 
N49 (e.g., Ichu fc Kennicutt 19881 : banas et al. 19971 : IrGose et al. 20041 : badenes et al. 20091 ). If 
N49 is a remnant of a core-collapse explosion from a massive progenitor rather than a Type la, the 
Si-rich nature of the Head and East regions may be generally considered to be explosive O-burning 
or incomplete Si-burning products from deep inside of the core-collapse SN. I n fact, Si-rich ejecta 
knots have been detected in core-c ollapse SNRs: e.g., Shrapnel A in Vela SNR (jMiyata et al. 200 ll ) 
and those found in Cassiopeia A (jHughes et al. 2 000). 



While the possibility of Type la origin for N49 is intriguing, we caution that the utility of our 
Mg;/Ms estimate to identify the SN type is limited, because it is based on small localized ejecta 
features whereas SN nucleosynthesis model calculations are for the integrated ejecta material from 
the entire SN. Thus, based on the current data, our discussion on the origin of N49 is far from 
conclusive. A more extensive ejecta search and comprehensive nucleosynthesis study are required 
to reveal the true origin of N49. We point out that the correct identification of N49's origin 
is particularly important because of its astrophysical implications in the context of the SNR's 
environment. For instance, if N49 is the remnant of a Type la SN, the long-standing argument for 
its physical association with SGR 0526-66 is unambiguously ruled out. A Type la origin for N49 
may also suggest an intriguing case for a prompt population la SN from a relatively young white 
dwarf progenitor in an alogy to the scenario suggested for SNR 0104-72.3 in the Small Magellanic 
Cloud (|Lee et al. 201 lh . 



-14- 



4.2. SGR 0526 66 

Our deep Chandra observation firmly establishes the previously suggested two-component X- 
ray spectrum for SGR 0526-66. The PL index (T ~ 2.5) of its X-ray spectrum is intermediate 
between those for other SGRs (T ~ 1-2) and AXPs (r ~ 3-4) which suggests that SGR 0526-66 



may be a transition object between the SGR and AXP states (jKulkarni et al. 20031 ). In fact, our 



estimated PL index (r ~ 2.5) is identical to that of AXP IE 1048.1- 5937, which is considered 



to be a clear case of a SGR- AXP transition object (jGavriil et al. 20021 ). The implied size of the 
BB radiation (i?BB ~ 5-6 km) is smaller than the canonical size of neutron stars. The small Rbb 
suggests the existence of restricted hot spots on the surface of the neutron star. 

Alternatively, a BB+BB model can equally describe the observed spectrum of 0526-66. The 
two-component BB model was preferred for some AXPs, because the second BB component is 
physically more r easonable than the steep P L component to explain the low flux limit at longer 



wavelengths (e.g., lHalpern &: Gotthelf 20051 ). Recently, Nakagawa et al. (2009) showed that the X- 
ray spectrum of the quiescent emission from 0526-66 can be fitted by a two-component BB model. 
Our result is consistent with that by Nakagawa et al. (2009). The hard BB component indicates a 
small emission area (R ~ 1 km) with a hot temperature of kT ~ 1 keV. The soft BB component 
indicates an area corresponding to the entire surface of the neutron star (R ~ 10 km) with a lower 
temperature of kT ~ 0.4 keV. These overall characteristics a re reminiscent of the peculiar types o f 



neutron stars found at the center of several young SNRs (e.g., iPavlov et al. 200(3 ; iPark et al. 20091 ). 
but the estimated BB temperatures of 0526-66 are significantly higher than those estimated in 
others. 

The overall X-ray flux in the 0.5-10 keV band in 2009 is ~15% lower than it was in 2000-2001. 
We show a long-term X-ray light curve of SGR 0526-66 in Figure [5j In Figure [5j we plot the mean 
flux of 2000 and 2001 data (the middle data point) because the observed flux is indistinguishable 
between the two epochs. We added the fluxes estimated by the ROSATKRI data in this light curve. 
For the R OS A T fluxes , we converted the observed HRI count rate (1.51+0.13 x 10~ 2 counts s , 



Rothschild et al. 1994 ) into the 0.5-10 keV band energy flux using PIMMS. In this calculation, we 
assumed two cases of BB+PL and BB+BB models with the best-fit parameters listed in Tables [4] 
and [5j respectively. We also assumed the fractional contributions in the observed 0.1-2.4 keV 
band HRI count rate from each of the model components, based on the results summarized in 
Table H and [51 The calculated ROSAT HRI fluxes are /0.5-10 keV ~ 1-34 (BB+PL) and ~1.17 
(BB+BB) x 10~ 12 erg cm -2 s _1 . While the X-ray flux appears to have decreased by ~30% since 
1992 based on the BB+PL case, the overall flux decrease is less certain for the BB+BB case. 
If it is real, the suggested X-ray flux change for 0526-66 would not be surprising, because long- 
term yariabilities by a factor of up to ^10 in se veral years have been detected in some AXPs 



(e.g., iBavkal et al. 19961 ; lOosterbroek et al. 19981 ). Follow-up Chandra observations are essential 



to reveal the true nature of the long-term light curve of SGR 0526-66. 



We note that the foreground column to 0526-66 shows a significant discrepancy by a factor 
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of ~3 between values estimated by two different spectral modelings (BB+PL vs. BB+BB). The 
column estimated by the BB+BB model fit (-/V"h,lmc ~ 1-7 x 10 21 cm -2 ) is generally in agreement 
with those measured for N49 (Table [2]) , and particularly, it is fully consistent with the column 
toward the SGR_BG region (Ah,lmc ~ 1-6 x 10 21 cm -2 ). These consistent columns between N49 
and 0526-66 are supportive of the long-suggested physical association between them. On the other 
hand, our BB+PL model fit of 0526-66 spectrum shows a substantially larger foreground column 
for SGR 0526-66 (Ah.lmc ~ 5.4 x 10 21 cm -2 ) than that for SGR.BG region. If this is the case, 
the large difference in the foreground absorption between SNR N49 and SGR 0526-66 brings into 
question their physical association. If we assumed an average interstellar density of no ~ 1-2 cm~ 3 
(see Section EE} near N49, SGR 0526-66 may be ~500-1000 pc beyond N49. Thus, this model- 
dependent difference in Ah,lmc for 0526-66 deserves full attention for further studies with follow-up 
observations. An H I survey of the LMC sho ws that N49 and 0526-66 are projected toward the 
dense boundary between two supergiant shells (|Kim et al. 20031 ). The estimated H I column density 
toward this region appears to be -/Vh,lmc ~ 5 x 10 21 cm~ 2 . It is difficult to discriminate our model- 
dependent A^h,lmc measurements toward 0526-66 based on this relatively large column estimated 
by H I data with a poor angular resolution (~1' which is comparable with the angular size of 
N49). We note that, if 0526-66 and N49 are associated, the projected angular separation (~22") of 
0526-66 from the geometric center of N49 requires a large kick- velocity (i.e., an average transverse 
velocity of v ~ 1100 km s _1 for the SNR age of ~4800 yr). At least, the proper measurement of 
the foreground column toward 0526-66 is directly related with two important astrophysical issues: 
(1) the origin of the quiescent X-ray emission from SGR 0526-66, and (2) the physical relationship 
between SNR N49 and SGR 0526-66. It may also provide a useful piece of puzzle to reveal the 
origin of SNR N49 (thermonuclear vs. core-collapse). 



5. SUMMARY 

Using our deep Chandra observation, we detect metal-rich ejecta features in SNR N49. These 
ejecta features include a "bullet" that is most likely a stellar fragment travelling beyond the south- 
western boundary of the SNR. We also find an ejecta feature in the eastern part of the SNR, 
nearly in the opposite side of the SNR to the bullet. Both of these ejecta features show highly 
enhanced Si and S abundances by an order of magnitude above the LMC values. We do not find 
compelling evidence for overabundant O and/or Fe in these ejecta features. If N49 is a remnant of 
a core-collapse explosion of a massive progenitor, the Si-rich nature of ejecta may be considered to 
be nucleosynthesis products of explosive O-burning or incomplete Si-burning from deep inside of 
the SN. On the other hand, the estimated Si- and S-rich ejecta mass ratio appears to favor a Type 
la origin for N49. If N49 was the remnant of a Type la SN, the suggested physical association 
between N49 and SGR 0526-66 would be unambiguously ruled out. However, we note that our 
SN ejecta study is limited because we used only some small localized ejecta features rather than 
the integrated ejecta composition from the entire SNR to represent the true SN nucleosynthesis. 
Follow-up studies to reveal the comprehensive ejecta composition in N49 are required to unveil the 



-16- 



true origin of this SNR. 

Our new Chandra observation allows us to detect the blast wave forming the swept-up shell in 
the southern boundary with significant photon statistics. Since this part of the shell is not affected 
by spectral complications due to the SNR's interaction with dense clumpy clouds, it provides a 
useful opportunity to reveal the dynamics of the SNR more accurately than previous works. We 
estimate the Sedov age of rs e d ~ 4800 yr and the explosion energy Eq ~ 1.8 x 10 51 erg for N49. 

Our spectral analysis of the quiescent X-ray emission from SGR 0526-66 using the deep expo- 
sure clearly reveals the presence of a BB emission (/cTbb ~ 0.44 keV) in addition to a PL component. 
The implied BB emitting area is relatively small (R ~ 5-6 km) compared to the canonical size of 
neutron stars. The estimated PL photon index (r ~ 2.5) is identical to that of AXP IE 1048.1- 
5937, the well-known candidate transition object between AXPs and SGRs. Alternatively, the 
observed X-ray spectrum of 0526-66 can be equally fitted by a two-component BB model. This 
model indicates that X-ray emission originates from a small (R ~ 1 km), hot (kT ~ 1 keV) spot(s) 
in addition to a cooler (kT ~ 0.4 keV) surface of the neutron star (R ~ 10 km). 

We find marginal evidence for a slow decay in the observed X-ray flux of 0526-66 (~20-30% 
for the last ~17 yr). Continuous X-ray monitoring of 0526-66 is needed to clarify the nature of 
its long-term light curve. We find a considerable difference in the foreground column by a factor 
of ~3 between two modelings (BB+PL vs. BB+BB) of the X-ray spectrum of 0526-66. While 
such a model-dependent discrepancy has been noticed for several other AXPs, the 0526-66 case 
is particularly intriguing because discriminating these competing models may be able to provide 
some critical clues on the nature of 0526-66 and N49: e.g., their physical association, the origin of 
X-ray emission of 0526-66, and the origin of N49. 

This work has been supported in part by the SAO under the Chandra grants GO9-0072A to 
the University of Texas at Arlington, and GO9-0072X to the Pennsylvania State University. J.P.H. 
has been supported in part by the Chandra grant GO0-11090X. 
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Table 1. Chandra Observations of N49 



ObsID 


Observation Date 


Exposure (ks) 


Instrument 


10123 


2009-7-18 


26.8 


ACIS-S3 (1/4 subarray) 


10806 


2009-9-19 


26.5 


ACIS-S3 (1/4 subarray) 


10807 


2009-9-16 


26.0 


ACIS-S3 (1/4 subarray) 


10808 


2009-7-31 


28.7 


ACIS-S3 (1/4 subarray) 


747 


2000-1-4 


39.9 


ACIS-S3 (1/8 subarray) 


1957 


2001-8-31 


48.4 


ACIS-S3 (1/8 subarray) 



Note. — In the Chandra archive, there are two other ObsIDs (1041 and 
2515) that detected N49. Those observations are not useful for the purposes 
of this work because of a large off-axis pointing (~6"5 for ObsID 1041) or 
a short exposure (7 ks for ObsID 2515). Thus, we excluded them in this 
work. 
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Table 2. Summary of Spectral Model Fits to Subregions of N49 



Region 


■Wh,lmc 
(10 21 cm" 2 ) 


kT 

(keV) 


n e t 

(10 11 cm" 3 s) 


(10 57 cm" 3 ) 


X 2 / n 


Shell 
SGR.BG 
Tail b 
Tail c 
Head 
East d 


0.89+^ 

-L-OO-0.44 
i ' ui -1.00 

i 09+0.60 

L - OZ -0.57 
1 1+O.8I 
°- li -0.88 


n r; 7 +0.05 
u -°' -0.10 

0.56±0.03 

2 n2 +1 - 41 

z - uz -0 .24 

n 71 + - 11 

1 nA+O-oe 

1 - u *-0.05 
1 nq+0.05 
i - uy -0.04 


D - oo -3.90 

9.65+H 

4 D9+ 3 - 78 
>20.3 
>67.0 


U - OU -0.28 

1 61 +0 - 40 
1 - Di -0.28 

9 9C-+0.19 

z -^ d -0.20 

2.68+ , : 2 ° 


49.1/68 
258.9/218 
44.9/44 
52.4/49 
67.0/59 
144.4/133 



Note. — Uncertainties are with a 90% confidence level. The 90% limit is presented 
where the bes t-fit value is unconstraine d. The Galactic column iVH,Gai is fixed at 0.6 
x 10 21 cm~ 2 (|Dickev fc Lockman 199ol ). 

a Volume emission measure, EM = J n e riHdV, assuming the distance to the LMC, 
d = 50 kpc. 

b The best-fit parameters are from a single shock model fit (see text). 

c The best-fit parameters are for the soft component from a two-temperature model 
fit (see text). 

d The best-fit parameters for the ejecta component. 



-21 - 



Table 3. Summary of Estimated Metal Abundances in Subregions of N49 



Region 





Ne 


Mg 


Si 


S 


Fe 


Shell 
SGR.BG 
Tail 
Head a 
East b 


n qo+0.49 
u -°^-0.18 

0.28^ 
0.32 

0.27 


27+ - 21 
u - z '-o.n 

17+ - 07 

u - i ' -0.06 

<1.54 (0.58) 
0.27 
0.21 


0.25±0.10 
1 5 +0 - 04 

a Q7+0.56 
-0.19 

0.25 
0.20 


3fi +0 - 17 

24+ ' 06 
u ' z *-0.06 

°- 88 -0:38 
9 98 +0.48 
Z - ZO -0.39 

1 Q2+ 029 


0.841°;! 
0.59l°;l 
<2.30 (0.49) 

o 9 +0.90 
°- z ^-0.80 

1 QO+0.44 
± - oy -0.39 


21+ - 12 
u - zl -0.08 

0.131°;°? 
0.2S$* 

0.21 

0.17 



Note. - Abundances are with respect to Solar (I Anders fe Grevesse 19891 ). Uncer- 
tainties are with a 90% confidence level. The 90% limit is presented where the best-fit 
value is unconstrained. 

a Abundances for O, Ne, Mg, and Fe are fixed at values measured for the "Shell" region. 

b Abundances for O, Ne, Mg, and Fe are fixed at the average values measured for the 
"Shell" and "SGR_BG" regions. 
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Table 4. Summary of PL+BB spectral model fit to SGR 0526-66 



Parameter BB PL Overall 

kT (keV) 

r 

i?BB,i a ( km ) 
#BB,2 b (km) 

/x,i a (10~ 12 erg cm' 2 s -1 ) 
fx,2 b (W~ 12 erg cm" 2 s" 1 ) 
L x ,i a (10 35 erg s- 1 ) 
L x , 2 b (10 35 erg s' 1 ) 
iVH.LMC (10 21 cm" 2 ) 
A^H.Gal (10 21 cm' 2 ) 
X 2 / n 



Note. — Uncertainties are with a 90% confidence level. 

a This is based on the Chandra data taken in 2000 and 2001. fx 
and Lx are estimated in the 0.5-10 keV band. The distance to the 
LMC, d = 50 kpc, is assumed. 

b This is based on the Chandra data taken in 2009. fx and Lx 
are estimated in the 0.5-10 keV band. The distance to the LMC, d 
= 50 kpc, is assumed. 



0.44±0.02 

2 50+ - 11 
z.ou_ 12 

5.5±0.6 

0.40±g;^ 0.78±0.13 1.18±g;lg 

0.33±g;gf 0.69±0.11 1.01±g;?| 

l-68±g;^ 3.84±0.65 5.52±g;^ 

1.3T±g"|| 3.36±0.54 4.74±g"|f 

0.44_ Q59 

0.6 (fixed) 
577.4/582 
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Table 5. Summary of BB+BB spectral model fit to SGR 0526-66 



Parameter BB so f t BBhard Overall 



kT (keV) 


0.39±0.01 


1 m+ 011 

i - ui -0.09 




i?BB,i a (km) 


q 7+0.6 


1.0±0.2 




i?BB,2 b (km) 


9.0±0.5 


0.9±0.2 




/x,i a (10~ 12 erg cm~ 2 s -1 ) 
/x,2 b (10 -12 erg cm" 2 s" 1 ) 
L x ,i a (10 35 erg s- 1 ) 
L x , 2 b (10 35 erg s" 1 ) 
Wh,lmc (10 21 cm" 2 ) 


n yo+0.09 
u - 10 -0.07 

62+ - 34 

u - DZ -0.23 

90+ ' 11 

u - yu -0.09 
n 77+O.43 
u -' '-0.28 


40+ - 22 
U - 4U -0.15 

S5+ ' 19 

zlo+0.23 
U - 4Z -0.16 

S6+ - 20 


i iq+0.01 
1 - lo -0.05 

96+ - 01 

u - yD -0.05 
o q«+0.05 
°- yu -0.18 

oq+0.04 

1 70+ - 25 

L - ' u -0.23 


iVH.Gal (10 21 cm~ 2 ) 






0.6 (fixed) 


X 2 / n 






616.4/582 



Note. — Uncertainties are with a 90% confidence level. 

a This is based on the Chandra data taken in 2000 and 2001. fx 
and Lx are estimated in the 0.5-10 keV band. The distance to 
the LMC, d = 50 kpc, is assumed. 

b This is based on the Chandra data taken in 2009. fx and Lx 
are estimated in the 0.5-10 keV band. The distance to the LMC, 
d = 50 kpc, is assumed. 
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Fig. 1. — (a) A 3-color ACIS images of N49 based on our 2009 data. Regions used for our spectral 
analysis are marked, (b) A zoom-in image of the SGR 0526-66 region. The source and background 
regions used for the spectral analysis of SGR 0526-66 are shown, (c) A zoom-in of the southwestern 
region of N49 including the ejecta bullet. Contours for the outer boundary of N49 are overlaid. In 
(a), (b), and (c), Color codes are: red is 0.3-0.8 keV, green is 0.8-1.7 keV, and blue is 1.7-7.0 keV 
band. The pixel size is 0"5 in all panels. 



5:26:000 550 25:50.0 10.0 05.0 5:26:00.0 55.0 

RA (J2000) RA (J2000) 



Fig. 2. — (a) An X-ray 3-color image of N49 with contours of Si EW map overlaid. Color codes are 
the same as those in FigureHJ The Si Hea + Lya (E = 1.75 - 2.1 keV) EW has been calculated by 
the methods described in Park et al. (2003a). (b) The same X-ray color image as in (a), overlaid 
with contours of an archival Hubble Space Telescope image. 




-5 - 



0.5 12 0.5 1 2 

Energy (keV) Energy (keV) 

Fig. 3. — (a) X-ray spectrum of "Shell" region, (b) X-ray spectrum of "SGR_BG" region. In (a) 
and (b), the best-fit plane-shock model is overlaid. The lower panel is the residuals from the best-fit 
model. 
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Fig. 4.— X-ray spectrum of (a) "Head", (b) "East", (c) "Tail", and (d) SGR 0526-66. In (a), (b), 
and (c), the best-fit plane-shock model is overlaid (solid curve). In (d), the best-fit BB+PL model 
is overlaid (solid curve). In (b), dashed and dotted curves show model components for the ejecta 
and the superposed emission from the shocked cloud, respectively. In (d), dashed and dotted curves 
show the best-fit PL and BB model components, respectively. The lower panel is the residuals from 
the best-fit model. 
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1.6 



1.4 



1.2 



0.8 



■ BB + PL model 
□ BB + BB model 



□ 



□ 



4.8x10" 



5xl0 4 



5.2xl0 4 
MJD 



5.4xl0 4 



5.6xl0 4 



Fig. 5. — The long-term X-ray light curve of SGR 0526-66. Fluxes at the first epoch have been 
estimated by the ROSAT HRI data ([Rothschild et al. 19941 ) . la error bars based on the count 
statistics are shown in the ROSAT fluxes. The next two fluxes are estimated by Chandra data. 
Two model fits (BB+PL and BB+BB) are assumed for the flux estimates. For the Chandra fluxes, 
uncertainties are with a 90% confidence level as estimated by "flux error" command in the XSPEC 
based on the two-component model fits. 



